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Abstract. In this paper we study the tailoring of photon spectral properties 
generated by four- wave mixing in a birefringent photonic crystal fibre (PCF). The 
aim is to produce intrinsically narrow-band photons and hence to achieve high non- 
classical interference visibility and generate high fidelity entanglement without any 
requirement for spectral filtering, leading to high effective detection efficiencies. We 
show unfiltered Hong-Ou-Mandel interference visibilities of 77% between photons from 
the same PCF, and 80% between separate sources. We compare results from modelling 
the PCF to these experiments and analyse photon purities. 
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1. Introduction 

Single photons are at the heart of optical quantum technologies such as quantum 
cryptography [I], quantum computing [2 J and quantum metrology [3|. These 
applications mainly rely on non-classical interference, so one constraining factor is that 
the photons must be in a pure quantum state [4]. A way to generate single photons is 
to use pair photon sources based on parametric down- conversion or four-wave mixing 
(FWM), where the detection of one photon of the pair can herald the presence of the 
second photon [5j. However, due to energy and momentum correlations, the heralded 
photon usually ends up in a mixed state [6]. The common way to counteract this issue 
is then to strongly filter the photon pair at the cost of losses, which significantly reduces 
the overall efficiency fi of detecting a photon. In a multi-photon experiment involving 
n photons (such as a cascaded CNOT gate), the overall probability of detection scales 
as [i n meaning that a high effective detection efficiency is essential. It has recently been 
suggested that naturally narrow band phase matching can be achieved in the case of 
a birefringent photonic crystal fibre (PCF) when the photon pairs are created in the 
orthogonal polarization to that of the pump light [7] . In this paper we study the tailoring 
of photon spectral properties generated by four-wave mixing in such birefringent PCF. 
The aim is to produce intrinsically narrow-band photons and hence to achieve high 
non-classical interference visibility without any requirement for spectral filtering. We 
examine experimentally how the phase-matching of several different PCFs compares to 
theoretical expectations, and perform non-classical interference to measure the purity- 
first between photons from two separate PCFs, and then from a single PCF pumped in 
two directions using a Sagnac loop configuration. We then use simulations to look at 
how the purity of the photons' state depends upon the fibre parameters and the pump 
laser and compare these to experimental results. Finally, we use the Sagnac setup to 
demonstrate a PCF source of polarization-entangled photon pairs, and perform quantum 
state tomography on the generated photons to quantify the entanglement present. 

2. Four- Wave Mixing In Birefringent Fibres 

At the quantum level the FWM process can be regarded as the virtual absorption of 
two pump (p) photons and subsequent creation of a signal (s) and idler (i) photon pair. 
The output state of the pair is then given by 



where the joint spectral amplitude (JSA) F(u s , uji) is equal to the product of the pump 
amplitude a(u s + ui) and a phase matching function <p(u s ,LUi). 

The process then has to satisfy energy conservation 




(1) 



Auj = 2lo p — uj s — uii = 



(2) 
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Figure 1. Signal and Idler wavelength as a function of Pump in a Birefringent fibre 



and a significant build up of amplitude in a length of fiber requires phase matching 

n v uj v n s oo s riiLOi n „ „ 
Ak = — - - 2 7 P P = (3) 



where n PtSy i are the refractive indices of the medium at pump, signal and idler 
wavelengths, P p is the peak pump power, and 7 is the nonlinear coefficient of the fibre: 



7 = 27m 2 /X p A eff (4) 

eff 



with n 2 = 2 x 10 20 m 2 /W the nonlinear refractive index of silica, and A e ft the effective 



cross-sectional area of the fibre mode [Sj. 

For a birefringent fibre, the refractive index will be different along the slow (s) and 
fast (f) axes. Solving equations[2]and[2]will result in four different solutions as illustrated 
in Fig. [TJ Either the pair is created along the same axis as the pump photons (ss — > ss 
and // — » //) or the pair is created along the axis orthogonal to that of the pump 
(ss — > ff and // — > ss). In the fabrication process of the fibre, a careful choice of 
n p and rij can lead to a situation where for an appropriate X p the group velocity of 
the pump is equal to that of the idler [6]. This situation corresponds to the zero slope 
region of the ss — » // curve in Fig. [1] giving rise to signal photons that are intrinsically 
narrowband and insensitive to pump wavelength over a short range. In this case the 
photon-pairs are created in an uncorrelated spectral state with a JSA that is factorable: 
F(u s ,Ui) = S(u) s ).I(cui), and the detection of one photon of the pair can then be used 
to herald the other in a pure quantum state. 
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3. Experimental study 

We have fabricated a set of birefringent PCFs (-A, -B, and -C) designed to have a zero 
dispersion wavelength (ZDW) around 780 — 800 nm. These fibres have similar dispersion 
profiles, but exhibit a slightly different birefringence due to different air hole size in the 
cladding as can be seen if Fig. [2j The ZDW and the birefringence were chosen to satisfy 
the cross polarization phase matching ss — > ff in the pumping region 700 — 740 nm 
accessible by our Ti: Sapphire laser. The zero slope section of the signal phase matching 
shown previously in Fig. [I] was made to sit in the region 600 — 650 nm. In this manner, 
the matching idler photons are emitted in a range where Si APDs are still efficient 
(800 — 870 nm), but are far enough from the pump wavelength to keep the Raman 
background low. 




Figure 2. Scanning electron microscope (SEM) image of the central region of PCF A. 
The birefringence required was introduced by reducing the size of two holes adjacent 
to the core (11 o'clock and 5 o'clock). 

Using a plane- wave method [9] , we can calculate the dispersion curves and retrieve 
the phase matching solutions for the three fabricated PCFs. We show in Fig. [3] 
the theoretical phase matching for the ss — > ff cross polarization case and the 
corresponding experimental results when pumping from 700 nm to 740 nm, which are 
in good agreement. 

3.1. Spectral bandwidths 

We now focus on the zero slope region of the phase matching, and study the spectral 
bandwidth of the signal photons as a function of fibre length. As described in [TO] , the 



Intrinsically narrowband pair photon generation in micro structured fibres 



5 




i -| — i — i — i — i — i — i — i — i — i — i — . — i — . — i — . — i — . — i — i — | 

670 680 690 700 710 720 730 740 750 760 770 

Pump wavelength (nm) 



Figure 3. Theoretical phase matching curves and the corresponding experimental 
results in the ss — > ff cross polarization phase matching case, for a set of three 
birefringent PCFs (-A, -B, and -C). 
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Figure 4. Spectral bandwidth of the signal photons as a function of fibre length and 
a 1/L fitting curve. 



bandwidth equation of the generated photons is given by 



Aw., 



2vrc 



\K-M\L 



+ 2 



M s - Mi 



(5) 



where L is the fibre length. 
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This equation is thus driven by two terms, one corresponding to the natural 
bandwidth of the phase matching, and a second related to the pump bandwidth 
broadening. The first is proportional to the inverse of the fibre length while the second 
is proportional to the slope of the phase matching curve. In our case, sitting on the zero 
slope region of the phase matching curve of the signal photons means that the signal 
bandwidth is insensitive to pump bandwidth, and thus only behaves as the inverse of 
the fibre length. This is verified by our measurement shown in Fig. H] where we can see 
the experimental data and the corresponding 1/L fitting curve. This study was limited 
by the fact that homogeneity over the length of fibre is very difficult to control during 
the fabrication process and lm of fibre is the longest usable length. As we will see in 
the next section, this limitation will also restrain our experimental study on purity. 

3.2. Quantum interference 

In order to experimentally test the purity of our heralded single photons, we performed 
a Hong-Ou-Mandel (HOM) experiment [11] using heralded single photons generated 
from independent PCF sources. Fig. [5] shows the setup of the experiment. Amongst 
the three PCFs fabricated, we decided to use PCF-A which has the highest birefrin- 
gence. In this manner, pumping at 705 nm corresponding to the zero slope region of the 
signal curve, phase matching is satisfied for generation of signal and idler photons at 
597 nm and 860 nm respectively. At these wavelengths, the detection efficiency is high 
(~ 60% for the signal and ~ 40% for the idler) and the Raman background in the idler 
arm remains low. We use a mode-locked TkSapphire pump laser (TiSa) emitting 100 fs 
pulses at 705 nm with a repetition rate of 80 Mhz which passes through a tuneable filter 
allowing control of the bandwidth from 7nm (~ 100 fs) to 0.25 nm (~ 3ps). Within 
this full range, the signal bandwidth remains constant at ~ 0.15 nm demonstrating its 
independence to pump wavelength in this region. The light is then sent onto a 50:50 
beam-splitter (BS) to pump two separate 40 cm long pieces of PCF-A. At the output of 




Figure 5. Schematic of the experimental setup to demonstrate quantum interference 
between photons generated in separate fibres. 
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Figure 6. HOM interference visibility plotted as a function of pump bandwidth. 

each PCF, a polarizer (Pol) oriented with the polarization of the pairs is used to filter 
the remaining pump light in the orthogonal polarization, and a dichroic mirror (DM) 
centered at 700 nm is used to split signal and idler photons into two different arms. 
In order to get rid of the residual pump light, two bandpass filters (BF) of 40 nm and 
10 nm bandwidth are used in the signal and idler arms respectively. The idler photons 
of each source are then launched into single mode fibres and the signal photons into a 
single-mode 50:50 coupler with polarization controllers (PC) in one of the input arms. 
As the photons are well mode matched to these optical fibres and they have not under- 
gone narrow spectral filtering we are able to achieve high overall lumped efficiencies for 
both sources of r] s = 0.25 and rji = 0.18 for signal and idler respectively. All four out- 
puts are finally connected to Silicon avalanche photodiodes (APD) linked to a four-fold 
coincidence measurement apparatus (&). 

A retro-reflector mounted on a micrometric translation stage in one arm of the 
interferometer allows us to vary the temporal overlap (5t) between the two photon 
pairs. When the two signal photons are made to arrive simultaneously on the beam 
splitter (St = 0), we can see a reduction of the four-fold coincidence count rate, which 
gives a measure of the purity of our heralded state. Using the tuneable filter, we 
have performed this experiment for different pump bandwidths. The visibility of the 
HOM dip as a function of pump bandwidth is plotted in Fig. |6] The visibility seems 
reasonably constant across the range shown with a weak peak when the pump bandwidth 
is ~ 3.03 nm where a 80 ± 2.5% visibility dip can be observed (see Fig. [7]). Although a 
clearly non-classical result, the visibility falls below the theoretical maximum of 100%. A 



Intrinsically narrowband pair photon generation in micro structured fibres 



8 




Figure 7. Coincidences as a function of delay position with a pump bandwidth of 
3.03 nm. 



possible reason for this is that the separate fibres produce signal photons with slightly 
different spectra due to variation in the fabrication process. To investigate this, we 
performed an experiment where both photons are created in the same fibre, and are 
thus spectrally indistinguishable. 



3.3. Quantum interference from a bi-directional source 

To test whether this limit to the visibility of the two photon interference was due to 
the use of two separate lengths of PCF we decided to try a method of generating two 
pairs of photons in one length of fibre, as was previously carried out in reference [12J. 
In this scheme the pump pulse is separated into two pulses that are set to counter- 
propagate through opposite inputs of the same fibre, but on the same principle axis. 
This is achieved through a Sagnac loop configuration which is also useful in generating 
entangled photons. 

The setup for this experiment can be seen in Fig. [HJ A 20 cm length of the 
birefringent PCF used previously is set up on two alignment stages with the ends able 
to rotate about the axis of propagation, allowing the fibre to be aligned along either 
birefringent axis. A polarizing beam splitter (PBS) is then placed at the point where 
the output modes of both ends of the PCF meet. It can be seen at this point that if the 
pump pulse is sent through the PBS it can travel to either input of the PCF dependent 
upon it's polarization. The PCF input for pumping clockwise around the loop is rotated 
so that the slow axis is horizontal, while the input for pumping counter clockwise has 
the slow-axis orientated vertically. 
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Figure 8. The setup used to generate two pairs of photons from one PCF and perform 
a polarization based HOM dip between a photon of each pair. 



Pump photons propagating in a counter clockwise direction are horizontally 
polarized and therefore these photons generate vertically polarized signal and idler 
photons due to the cross-polar phase matching scheme. As the fibre is twisted by 
90° the signal and idler polarizations are rotated to be horizontally polarized and exit 
through the PBS towards the detectors, leaving the output state as \^)) ou t = \H) s \H)i. 
If the pump photons propagate in the clockwise direction then the same processes occur 
but H — > V and V — > H. The vertically polarized pump photons \ip)i n = \V)p 1 \V)p 2 
generate horizontally polarized signal and idler photons \ip) = \H) s \H)i which are then 
transformed back to vertical polarization by the twist in the fibre so that \ip)out = 
\V) s \V)i. These photons therefore exit again through the same output of the PBS as in 
the counter clockwise case. Any pump photons which do not generate signal and idler 
photons are rejected back towards the pump laser. 

If the input power is large enough there is a finite probability of generating a pair 
of photons in both directions, so we have the state 

l^out = \H) 9 \H)i\V) 9 \V)i (6) 

all in the same single spatial mode. By inserting a dichroic mirror we can split the signal 
and idler photons into separate modes for heralding and analysis. All of the photons 
are collected into single mode fibres with polarization controllers to correct for any 
rotation caused by bends in the fibre and sent to analysis benches. The idler photons 
are split into two polarization modes by a PBS and then detected using silicon APD's 
as before, heralding the signal photons in the setup. The signal photons are passed 
through a HWP in a rotation mount and then sent through another PBS, followed by 
detection. Here we use the HWP in the signal arm to perform a polarization analogue 
of the standard delay-based HOM dip, where the polarizations of the signal photons are 
rotated from H and V to +45° and —45° respectively, and the degree of suppression of 
the heralded (and therefore four-fold) coincidence detection rate across all output ports 
characterizes the purity and distinguishability of the incident photons [12]. With a pump 
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Figure 9. Interference fringes in the coincidences detected as the polarization is varied 
with a HWP. 

bandwidth of ~ 0.7 nm and signal and idler filtering the same as for the previous HOM 
dip experiment, we see in Fig. [9] fringes of ~ 77% visibility, again below the theoretical 
maximum of ~ 100%. To discern whether there is an inherent purity limitation with 
the PCF used we construct a model to compare to experiment. 

4. Schmidt number optimization 

The purity of the heralded state is determined by the degree of factorability of the 
JSA [13] and varies with pump wavelength, fibre length, pump bandwidth and fibre 
geometry (by which the group velocity for the signal, idler and pump wavelengths can 
be controlled). We have developed a piece of software to optimise these variables for 
maximum purity. A general JSA can be decomposed into its Schmidt modes 

F(u„Ui) = ^2 \^jfj( u s)9j(ui) (7) 

3 

where fj(u s )gj(ui) are an orthonormal set of completely factorable modes. To quantify 
the factorability, we use the singular value decomposition (SVD), which is the matrix 
analogue of the Schmidt decomposition. We define the matrix representation of the 
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JSA as M = A$, where A is the matrix corresponding to the pump amplitude and 
$ to the phase matching function. The SVD is then a factorisation of M into three 
matrices, which will be called I,S,D: the columns of I represent the part of the Schmidt 
modes relating to the idler while the rows of S represent the modes of the signal, and D 
is a diagonal matrix linking the two. The elements of the diagonal matrix contain the 
amplitudes and are in descending order of magnitude. 

We obtain D from the SVD and calculate the Schmidt magnitudes Xj, with the 
following normalisation: 

5> = 1 (8) 

3 

A completely pure state will be represented by \j=i = 1 and A^i = 0. The Schmidt 
number K is then defined as 

For a pure state, K will equal 1. For a completely entangled state, there will 
be an infinite number of Schmidt magnitudes A, each approaching zero. The sum of 
the magnitudes will therefore approach zero with the result that K — » oo. We have 
previously seen that the JSA is defined as F(cu s ,cOi) = a(u s + Ui)(p(Lo s ,LOi), where 

a(u s + u)i) = e * (10) 

and 



iKL 



,AkL, 



e 2 sinci - ) (11 



Following the same steps described in [TO], we develop the phase- matching equation 
around the phasematched frequencies using a first order Taylor series in order to express 
A A; in terms of the group index. 



Ak 



2k 
2 



p 

dk 



du 



k s k{ 



2 7 P 
dk s 
dui 



(12) 



Auu — 

s du 



Au 



With the help of energy conservation, we can substitute Au p and find: 

(hft-KTA 

Ak = ^ ^Aoo ks + 



(A/" s — M p ) (Afj-Afp) 

Ao;fc s H 

c c 



Aui 



(13) 



where A/j = uji ^ | ^ + ni is the group index for mode 1. The JSA, and therefore M can 
now be expressed as a function of Au s and Aui only. Using the Sellmeier equation p3] , 
we can then generate the group index for a given fibre geometry and study the factora- 
bility of the JSA for a given phase-matching function, pump bandwidth and fibre length. 
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Figure 10. (a) Schmidt number of the heralded state as a function of pump 
bandwidths for three different fibre lengths: 20, 40, and 80 cm. (b) Minimum Schmidt 
number and the corresponding pump bandwidth as a function of fibre length. 





Figure 11. JSA of the created photon pair given by the product of the phase matching 
function and the envelope of the pump intensity (a) Unfiltered case (b) Where large 
bandpass filters (AXfn ter 3> signal/ idler ) have been applied to both photons. 



We first look at the purity of the heralded state as a function of pump bandwidths 
for three different fibre lengths: 20, 40, and 80 cm. The result of the simulation is 
shown in Fig. 10(a) As we can see, each fibre length has an optimum pump bandwidth 
corresponding to a minimum of the Schmidt number. We then plot in Fig. 10(b) the 
minimum Schmidt number accessible and the corresponding optimum pump bandwidth 
as a function of fibre length. This shows that the longer the fibre the better purity (lower 
Schmidt number) we can get. It is also interesting to note that the optimum pump 
bandwidth decreases when the fibre length increases. One explanation is that when 
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the pump bandwidth is wide, both signal and idler are limited in width by the phase- 
matching function, while when it is narrow, both are limited by energy conservation. 
At the optimum bandwidth the signal is mainly limited by the phase-matching and 
the idler mainly by energy-conservation, so their state is more factorable. As the fibre 
length is increased, the phase- matching sine function decreases in width (equation ITTj) . 
and so the pump bandwidth can also be decreased without reaching the region where it 
limits the signal bandwidth. 




Figure 12. Unfiltered Joint Spectral Intensity for a length of AO cm and pump 
bandwidth 2.6 nm. 




Wavelength (nm) Wavelength (nm) 



(a) (b) 
Figure 13. Spectral intensities projected from the JSI for (a) the signal (b) the idler. 



Fig. 10(b) shows that even in the case where the pump bandwidth has been 
optimised for a given fibre length, the Schmidt number still does not reach the optimum 
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value of 1. Looking at the JSA for an example of fibre length of 40 cm (Fig. 11(a) ), we 
can see that there are still some correlations between the signal and idler frequencies 
coming from the residual ripples of the phase matching sine function. These ripples 
will appear small in the measurable spectrum, shown as a Joint Spectral Intensity (JSI) 
in Fig. [12] and spectral intensities of the signal and idler in Fig. [131 but the resulting 
correlations in the JSA will increase the Schmidt number. A way to enhance the purity 
of the heralded photon is then to remove these ripples using bandpass filters whose 
bandwidths can be much greater than the full-width half maximum (FWHM) of the 
generated photons. This will result in an improvement of the purity without involving 



significant losses in photon collection efficiencies. Fig. 11(b) shows the same fibre 



length as in Fig. 11(a) where bandpass filters were applied to both signal and idler 
photon spectra and the pump bandwidth has been optimized to this new configuration. 
As a result, we can see that most of the ripples have disappeared and the calculated 
Schmidt number has decreased from 1.04 without filters to 1.01 in the filtered case. It is 
important to note that the photons can still be considered as unfiltered since the filters 
used satisfy the following condition: AX fUter > A\ signa i /idler . 

We expect the Schmidt number to be related to the HOM dip visibility by 
V = 1/K. To determine whether the HOM dip visibilities seen experimentally were 
limited by correlations in the JSA we can compare the Schmidt numbers derived from 
the experimental results with those obtained from the simulation for a range of pump 
bandwidths and the experimental conditions used. We show in Fig. [TH that there is 
an approximately constant offset between the experimental and theoretical results for 
the two fibre case, with the experimental Schmidt numbers never reaching the expected 
low values of the simulation, suggesting there is another reason for the lack of photon 
purity. To test whether this difference arises from the use of two separate pieces of 
fibre we can use the visibility of the fringes from the bidirectional case to estimate a 
Schmidt number for this fibre length and pump bandwidth, and also compare this to 
the simulation. In doing so we see in Fig. [TH that there is no offset present and the 
theoretical Schmidt number lies within the error of the experimental point, showing that 
the lack of visibility could be entirely explained by correlations between signal and idler 
photons being present in the JSA. This would confirm that the extra loss of visibility 
in the two fibre experiment is explained by differences between the two lengths of fibre 
introducing some cross-source distinguishability. 

We can also note that if the bidirectional source was pumped with a larger pump 
bandwidth of ~ 4 nm we can expect a fringe visibility of over 95% and a Schmidt number 
below 1.05. Further experimental investigation is being carried out to confirm this result. 
We note also that the pump filtering method produces top hat pump envelope shapes 
which increase the Schmidt number from that which is achieved through a Gaussian 
pump envelope of similar bandwidth. 
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Figure 14. Comparison of experimental and theoretical Schmidt number for a variety 
of fibre lengths. Black squares show two fibre experiment points while red circles show 
the bidirectional results. 



5. Entangled photon source 

The Sagnac loop configuration used above has another intrinsic use when used at lower 
powers, namely the generation of entanglement. The use of entangled photons is central 
to many quantum communication and computation schemes [U El [15] and plays a vital 
role in the development of quantum repeaters [16] and quantum memories |17j . The 
Sagnac loop method has been well demonstrated previously for degenerate and non- 
degenerate photons, both in fibre [T8l [T9] and using nonlinear crystals [20] , however all 
implementations required filtering and therefore suffer from low counting rates. Here 
we show a demonstration of a source of entangled photons which are in an intrinsically 
pure state from a PCF in a Sagnac loop, therefore requiring no narrow spectral filtering. 

The setup for this experiment is very similar to the bi-directional interference 
experiment with extra wave plates for analysis and a Soleil-Babinet compensator added, 
as seen in Fig. [15j A Ti:Sa laser emits 100 fs laser pulses (~ 8nm broad) that are 
then filtered as before, using two diffraction gratings and a slit, down to around lnm. 
This is then passed through a HWP which allows us to set a linear superposition of 
horizontal and vertical polarization followed by a Soleil-Babinet compensator consisting 
of two birefringent crystal wedges that can slide across one another, thereby imparting 
a variable phase between horizontal and vertical components. With this combination 
any pump state \ip)in can be set for the input of the system. 
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Figure 15. The setup used to generate pairs of entangled photons from a PCF in a 
Sagnac loop. 



\^)in = a\i)H + ^\i)v (14) 

Here we use 7 to signify a coherent state of many photons, as are present in the 
pump pulses. This pulse then enters a PBS which transmits horizontally, and reflects 
vertically, polarized photons. When two horizontally polarized photons travel in the 
counterclockwise direction around the 20 cm loop of PCF they can generate vertically 
polarized signal and idler photons through FWM in the cross-polar regime. These 
are then converted back to horizontal polarization by a 90° twist in the fibre about 
the axis of propagation, meaning that \H)p 1 \H)p 2 — > \H) s \H)i. When two vertically 
polarized pump photons travel in the clockwise direction they generate horizontally 
polarized signal and idler photons that are then twisted back to vertical so we have 
|y)p 1 |y)p 2 — >■ |V) s |V)i. The twist in the fibre also ensures that photons propagating in 
both directions are launched onto the same axis of the PCF, so they see the same phase 
matching conditions. If we then take into account the a and weightings set by the 
input HWP we have the following at detection 

|VW = a\H) 9 \H)i + (3e 2 ^\V) s \V)i (15) 

As the pump is orthogonally polarized to the signal and idler photons it is 
automatically rejected back towards the laser, reducing background noise at the output 
port. We can see from equation [15] that if we set a = (3 = and set = then 
we generate the maximally entangled |$ + ) Bell state. It is worth noting that this 
is an inherently stable setup in that the output state is not affected by temperature 
fluctuations or vibrations as the relative phase between photons generated in different 
directions is automatically compensated [2T]. The output state does not depend upon 
the length of the fibre or the point in the fibre where the photons are generated, only on 
the setting of the input state, making this a very versatile and useful source for quantum 
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information applications. At the output of the PBS the signal and idler photons are 
split to two paths by a dichroic mirror and then collected into optical fibres, each with a 
polarization controller to compensate for any stress induced polarization rotations, and 
then sent to analysis benches containing a QWP, HWP and PBS. The photons are then 
launched into multimode fibre and sent to silicon APDs on the transmitted (H) output. 

5.1. A test of entanglement 

An initial test of entanglement is a measurement of coincidence fringes. The analysis 
state of one photon of the entangled pair is set on one of two non-orthogonal bases, 
typically the computational {H/V} and diagonal {D/A} bases. The other photon is 
rotated continuously through the basis states and the coincidence rate between the 
detection of both photons is recorded. A second fringe is then measured for the 
other (non-orthogonal) basis. The classical {H/V} basis measurement should yield 
a coincidence fringe with unit visibility in the absence of all background noise, whilst 
the {D/A} basis measurement should be flat, or have zero visibility, for no entanglement 
and have unit visibility if the photons are maximally entangled. If the visibility in this 
basis is greater than 71% then the state is entangled enough to violate a Bell inequality 
[22| [23] . making this the lower bound on proving that the photons are entangled. 




signs! 

Figure 16. Coincidence count fringes taken with the idler photon fixed at 0° in black 
and 22.5° in red, showing 98% and 90% visibility respectively. 

For this test of entanglement the QWPs can be removed from the analysis benches 
for the signal and idler photons as we are only concerned with linear polarizations. We 
choose the idler photon to have the fixed polarization and set the HWP in this arm to 
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be 0°, or \H), for the {H/V} basis and 22.5°, or \D), for the {D/A} basis. For these 
two positions we rotate the HWP in the signal analysis bench whilst recording single 
APD counts and the coincidence between the signal and idler counts. The result in Fig. 
[TBI shows a visibility of 98% in the {H/V} basis and 90% in the {D/A} basis, clearly 
showing a high level of entanglement. 

The level of entanglement could be reduced by reflections in the fibre [21], excess 
background noise from residual pump light, Raman photons, or most likely from an 
imbalanced collection of photons from one direction of propagation to the other. This 
is difficult to correct for as when changing the output of one end of the loop you are 
changing the input coupling, and vice versa, which are often countering each other 
in quality. This is, however, not a fundamental limit to the amount of entanglement 
possible with such a setup, which should theoretically produce maximally entangled 
pairs. 

5.2. Quantum state tomography 

A further and more in depth analysis of the entangled state can be carried out to 
reconstruct the density matrix of the photon pair, known as quantum state tomography 
or QST [24J. To perform QST one must take many copies of the state to be analysed and 
perform a tomographically complete set of measurements, where the set of measurement 
operators define a basis in the Hilbert space of the two photon state. In this case that 
complete set of measurements is 16 polarization settings for the HWP's and QWP's in 



Table 1. Wave plate settings for two photon state tomography [24] 
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Figure 17. Tomography of the generated entangled state showing a.) the 
reconstructed real and imaginary parts of the density matrix and b.) the real and 
imaginary parts of the ideal $ + state. The experimental result shows a fidelity of 84% 
to the |$ + ) state. 



the signal and idler analysis benches. It comprises combinations of analysing \H), \V), 
\D) = ±(\H) + \R) = ±(\H) - i\V)) and \L) = j-(\H) + i\V)) states as shown 
in Table 1. 

Taking the 16 counts obtained for the different measurement settings we can 
perform a linear tomographic reconstruction of the state. This does, however, sometimes 
recreate states that are non-physical as they have negative eigenvalues and the trace of 
the square of the density matrix becomes greater than the allowed value of 1 . This means 
the density matrix p is not positive semi-definite as required. To avoid this problem we 
employ a technique, described in [23] and [25], known as maximum likelihood estimation. 
This approach estimates the density matrix by first generating a formula for a physical 



Intrinsically narrowband pair photon generation in micro structured fibres 



20 



density matrix and then using a likelihood function which quantifies how closely the 
estimated state would generate coincidences compared with the experimental data. This 
process is then repeated and, using standard numerical optimization techniques, the 
optimum physical density matrix is found. This process can also be improved by the 
use of convex optimization techniques [26]. Upon performing this calculation we find 
the closest density matrix to be 



Pexp 



( 



\ 



0.446 
-0.033 — 0.021* 
0.002 - 0.055i 
0.382 - 0.042i 



-0.033 + 0.021i 

0.049 
0.029 + 0.008i 
-0.053 + 0.064z 



0.002 + 0.055i 
0.029 - 0.008i 

0.037 
0.029 + 0.065i 



0.382 + 0.042i \ 
-0.053 - 0.064i 
0.029 - 0.065i 
0.467 J 



(16) 



which can be seen graphically in Fig. [17] split into it's real and imaginary parts. Also in 
Fig. [17] are the real and imaginary parts of the ideal |$) + state, which when compared 
to our experimentally generated state shows a fidelity of 83.93% ± 0.08%. 

We can also quantify the amount of entanglement generated by looking at the 
Tangle, T = 0.5495 ± 0.0027 which is greater than 0.5 giving a clear sign that 
entanglement is generated [23]. The Concurrence C = 0.7413+0.018 and Linear Entropy 
C = 0.3306 + 0.0013 also show a degree of entanglement but are clearly not the optimum 
values of C = 1 and £ = 0. As mentioned before, this non-maximal entanglement is 
probably due to the \H) and \V) modes not being perfectly overlapped. Nevertheless, 
this amount of entanglement is already at a useful level, and may be improved with 
careful alignment. 



6. Applications and Future Work 

Many further experiments are possible with this source. The bi-directional source set- 
up, for instance, can be extended to 2 photon NOON-state metrology, to show phase 
measurements with higher precision than allowed by the standard quantum limit, at 
very low intensities. The idler photons are used to herald the signal photons, which are 
injected into a displaced Sagnac interferometer, so that when the phase in one arm of 
the interferometer is varied, the probability of a coincidence across both outputs shows 
super-resolution of fringes with a spacing of 312 nm: half the wavelength of the signal, 
and less than half that of the pump. This compares favourably to similar experiments 
using spontaneous parametric down- conversion sources, where the fringe spacing is equal 
to that of the pump [27] . 

By using one entangled source and one single direction correlated source, a 
teleportation scheme [28] between photons of different wavelengths can be demonstrated. 
Here, the signal photon from the correlated source heralds the presence of the idler, 
which can be polarization-encoded to form a qubit using wave plates. When this 
photon and the idler from the entangled source are mixed on a 50:50 beamsplitter 
and a coincidence is detected at the two outputs, the two photons have been projected 
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onto a Bell state, and the qubit encoded previously is teleported onto the state of the 
remaining signal photon from the entangled source following a polarization rotation of 
45°. Teleportation has applications in quantum communications, for instance quantum 
key distribution [TJ, and in linear optical quantum computing [2J, where the success 
probability of a non-deterministic gate can be improved to an arbitrary level with the 
addition of teleportation and a resource of entangled ancilla photons. 

Similarly using one entangled source and a correlated source, a scalable version 
of a quantum CNOT gate can be performed [29], between a control and target qubit 
encoded in photons of different wavelengths. The signal and idler photons from the 
correlated source are used as the target and control qubits respectively, and an input 
state is encoded in their polarization using wave plate rotations. The control photon is 
mixed on a polarizing beamsplitter (PBS 1) with the idler from the entangled pair, and 
the target is mixed on a separate polarizing beamsplitter rotated by 45° (PBS 2) with 
the signal from the entangled pair. A successful CNOT operation has taken place if a 
diagonally polarized photon is detected at one output of PBS 1 in coincidence with a 
horizontally polarized photon at one output of PBS 2. The final state control and target 
photons are left at the remaining two outputs. The fact that a successful operation is 
heralded by ancilla photons, rather than relying on detecting the control and target as 
in [301 ED E21 [33j [3U [351 [36], means this CNOT is scalable and could be used as part of 
a more complex operation. 

Another application is in the generation of larger entangled states such as cluster 
states. In the cluster-state, or one-way, model of quantum computation, a multipartite 
entangled state is used as a resource for running a quantum algorithm represented by a 
series of single-qubit measurements [37] . Using two entangled sources and a polarizing 
beam splitter acting as a fusion gate we can generate a four photon cluster state and 
demonstrate quantum gates in the one-way model. Fusion gates are a resource efficient 
way of building up cluster statespj)], and we can extend this to larger states by adding 
more entangled pairs and fusion gates, or by encoding more qubits in other degrees of 
freedom, such as path or time-bin encoding. 

In conclusion, we have demonstrated a fibre-based source of photon pairs which 
avoids two significant limitations on the efficiency of conventional downconversion 
sources: the flat phasematching curve means the photons are created in an intrinsi- 
cally pure state without the need for spectral filtering, and the photons are created in 
a single-mode, so we achieve good spatial modematching to fibre. We have also shown 
high visibility non-classical interference between photons from the same and separate 
sources, and generated entangled pairs with a good fidelity to a maximally entangled 
Bell state. This makes the source suitable for many quantum information applications 
and future experiments. 
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